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1 .   Introduc'slon. 

Poisson's  integral  equation 

C3 

n  -n       rfvi     x     a  6(t)  dt 

—CO 

may  be  regarded  as  the  one  which  arises  in  connection  with  the 
problem  of  finding  a  potential  function  ^(x,y)  which  satisfies 

■>b       +  ib         =  0  y  >  0 

*xx   • yy  J 

and  which  is  equal  to  f(x)  along  the  line  y  =  a  >  0.   If 
V(x,0+)  =  <f>(x);  and  if  #(x,y)  is  regular  at  infinity  then 

CO 

If  ^(x,a)  =  £{::)    the  function  <f)(t)  must  satisfy  Poisson's 
equation.   The  function  f(x)  is  called  the  Poisson  transform 
of  <|>(t)  and  we  will  write 

CO 

f  (x)  =  P  <|>  s  1   .  f       a  tltl  dt_ 


a"  "  ^    (t-x)2  +  a2   ' 

—  CO 

We  can  say  that  Poisson's  integral  equation  poses  the 
problem  of  finding  the  boundary  value  of  a  function  known  to  be 
harmonic  in  a  half-plane  when  the  boundary  is  inaccessible,  but 
measurements  can  be  taken  along  a  line  parallel  to  the  boundary. 
This  kind  of  problem  arises  In  various  physical  investigations. 
For  example,  the  solution  of  the  problem  is  of  considerable  value 
for  the  interpretations  of  observations  derived  from  certain 
methods  of  geological  prospecting  which  depend  on  a  survey  of  a 
potential  field  such  as  the  gravimetric  or  magnetic  field  of  the 


. 


■ 


2 
earth.   In  prospecting,  the  problem  is  to  locate  and  evaluate 
concentrations  of  matter  which  cause  anomalous  fields  and  hence 
the  quantity  a  may  also  be  unknown  at  the  start. 

Snow's  formula  for  the  solution  of  (1.1)  is  well  known; 
but  not  long  ago,  Pollard  [1]  showed  that  the  Poisson  transform 
obeys  the  lav; 

(i.2)         pa  pb  <>  =  pa+b  * 

and  he  used  this  to  establish  Snow's  inverstion  formula  under  more 

general  assumptions  than  those  used  previously.  The  purpose  of 

this  report  is  to  present  other  proofs,  both  new  and  old,  of 

Pollard's  result;  and  to  show  some  extensions  and  applications 

of  the  property  (1.2). 

Section  2  contains  an  exposition  and  simplification  of 

Pollard's  proof  besides  some  other  methods  for  the  solution  of 

(1.1).   In  Section  5  it  is  shown  that  if 

P   2ir 
P   =  t1-*  )   /       <HcO  ch 

then 

?sl  A2       "A1A2 
and  this  is  used  to  solve  Poisson' s  equation  for  the  circle, 
namely 

f(G)  =  P^  <i>   . 

Section  4  is  concerned  with  a  generalization  of  Poisson' s  equation, 
the  integral  equation  which  arises  when  we  wish  to  determine  the 
boundary  value  of  a  two-dimensional  harmonic  functions  from  its 


7 
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behavior  alon^  an  interior  closed  path.   The  results  of  Section  1 
are  used  to  solve  a  certain  blharmonic  problem  posed  in  Section  5- 
In  Section  6   various  inversion  formulas  derived  in  Section  2  are 
examined  from  the  practical  point  of  view.  Several  approximation 
procedures  are  suggested  and  Section  6  includes  a  derivation  of 
an  iteration  procedure  for  the  solution  of  Poisson's  equation. 
Finally,  Section  7  shows  how  some  of  the  foregoing  results  can  be 
extended  to  analogous  potential  problems  in  higher  dimensions. 

2 .   Some  Methods  for  the  Solution  of  Poisson's  Integral  Equation . 

We  assume  that  <|>(£)  is  lntegrable  on  each  finite  interval, 
and  that  each  of 

(2.D        fiii^;      r±u)dt       c>0 

t  -■     t 

C  -w 

exists.      Under   this   assumption,    the   integral 

00  00 

*    J      (t-x)2+a2       2lrl  l-x-ia       t-x+^a 

-CO         v  — oo 

exists,  if  a  >  0;  and  it  represents  a  function  of  x  which  is 
analytic  when  x  is  real.  Hence  a  solution  of  Poisson's  integral 
equation 

00 

(2.2)        f(x]  -  -   /     ...2^2  a  >  0 

(t-x)  +a 

—  00 

cannot  exist  if  the  prescribed  f(x)  fails  to  be  differentiable 
an  infinite  number  of  times  at  any  value  of  x  which  is  real. 


. 
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I.   Introduce 

CO  00 

r-,/  \   1   /  a  <b(t)  dt    1    [  r      1        1   t  it^\    jf 
P(2)  =  1?   /   ,/  X    2  =  2H   '  [tTiTIi  -  TTi+iil  ^^  dt 
(t-z)  +a 

-00     X  —  CO 

where  z  =  x+iy.  This  function  of  the  complex  variable  z  is  analytic 
everywhere  except  along  the  lines  z  =  t-ia  and  z  =  t+ia,  each  of 
which  is  parallel  to  the  real  axis  and  a  units  distant  from  it. 
Since  F(z)  coincides  with  f(x)  when  z  is  real,  the  integral  which 
defines  F(z)  gives  the  analytic  continuation  of  f(x)  into  the 
strip  -a<  Im  z  =  y  <  a.   In  other  words,  if  i(t)  is  to  exist,  the 
prescribed  f(x)  must  be  analytic  on  the  real  axis  and  it  must 
allow  an  analytic  continuation  of  itself  into  the  strip  -a<Im  z  <  a. 
If  f(x)  allows  such  a  continuation,  the  continuation  is  unique. 

We  proceed  under  the  assumption  that  f(x)  satisfies  the 
analyticity  conditions  which  are  required  for  the  existence  of 
<t>(t)  .   If  e  >  0  we  have 

lim  (f[x-!-i(a-e]+f[x-i(a-e)]} 
e-o-0 


-  am   f1   r  £  Ht)   dt   1  p    (2a-e)  6(t)  dt, 
t-x  -i-o-   T  ■ 

—CO  —  CO 


e_>0  v  (t-x)N-o'-   ^     (t-x)4£+(2a~e)' 

—CO  —CO     * 

It  is  known  that  for  almost  all  values  of  x 


CO 


Li    \   i4   1   ('  e  4(t)  dt 
<J>(x)  =  lim  -   /  i-i-^ — n    . 

e->0  T  /   (t-x)2+e^ 


—  CO 


Hence  for  almost  all  values  of  x 


<f>(x)    lim  [f[x+i(a-e)]  +  f[x-i(a-e)]}  -i  f  2a  ^^  % 


s->0  T  -   (t-xr+4a 

—  00 

The  integral  in  the  last  equation  is  equal  to  ^(x,2a) 
where  ^(x,y)  is  the  potential  function 


' 


- 


- 
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_OT  (t-x)  +y 

According  to  (2.2),  this  potential  function  at  y  =  a,  namely 
tp(x,a),    is  prescribed  to  be  equal  to  the  analytic  function  f  (x)  . 
Consequently,  the  Poisson  integral  formula  for  a  function  harmonic 
in  a  half-plane  shows  that  ^(x,y)  is  given  by 

(2.3)  -7(x,y)  =|  f"   (ya)/(t)  dt 

T  _M  (t-x)"  +  (y-a)2 

when  y  >  a.   The  analyticity  condition  on  f(x)  insures  the 
uniqueness  of  #(x,y) .   It  follows  from  (2.3)  that 

CO  00 

1       '     2a  (|>(t)  dt  _   ,   2  ,  _  1   ;r  a  f(t)  dt 
w   _-;   (t-x)2+4a2  *  _^  (t-x^+a2 

We  therefore  conclude  that  for  almost  all  x 

CO 

(2.4)  t(x)  =  lim  (f[x+i(a-e)]  +  f[x-i(a-e)]}- -J  '  a   f(tj>  % 

e->0  u  -  (t-x)'+a2 

—CO 

is  the  solution  of  Poisson' s  integral  equation 

CO 

,.,  x   1  r  a  <j>(t)  dt  -co  <  x  <  co 

1 1  x ;  =  —  j  ^ — ^ 

u  j  (t-x  -Ki  a  >  0 

-co  x 

in  which  f(x)  and  <|>(x)  satisfy  the  conditions  stated  above. 
The  formula  (2.4)  can  be  expressed  in  a  slightly  more  compact 
and  theoretically  more  useful  form  if  we  use  the  fact  that  f (x) 
must  allow  an  analytic  extension  of  itself  into  the  strip 
-a  <   z  <  a.   If  we  start  with 

CD 

4>(x)  =  lim  (f[x+i(a-2e)]  +  f[x-i{a-2e)]J  -|  ,f_a_fjt)dt 
g->0  "    -'   (t-xr+a* 

—  CO 

on  what  is  the  same  thing, 

CO 

4>(x)   =  11m     2imi  f[x+i(a-2s)]    -^    Im    f    1{V   *'c   , 

£_>0  Tf  -'  t~-~la 

—CO 


' 


an  application  of  the  Cauchy  integral  formula  gives 

co  <»+i(a-£) 

Lt     ^  n        flTn,r  (t)     dt  1-      f  f(t)     rlt  , 

4>(x)   =   lira  {-  Im/      fc.;:.1(a.2e)   -  ¥  Im  /  t  -  :c  -  i(a-2s  ) } 

£">0  -«+i'(a-e) 

—  CO 


1  T«,  r  f(t)  c' 

-7   In7    t-x-a 


CO 

dt 

ia 


co+i(a-e) 
*(x)  =  -lim    llm/      t,^l(f!2£)    • 


~">0     -oo-KL(a-e) 

A   change  in  the  variable  of  integration  namely,    t  =  z  +  i(a-s), 
yields 

CO 

£->0 

-co 

or  its  equivalent 

C3 

»„  r\  \t    \        -n   Tm  1  r      f[z-i(a-e)]  dz 

(2.6)        <j>(x)  =  lim  Im-   /    Lz  _  x  ,  -£    • 

£->0 


—co 


II.  H.  Pollard  [1]  deduced  the  inverstion  formula  (2.4) 
for  the  Poisson  transform 

CO 

(2.7)       f<x)  - ±  ra*'Hdi  *  pa  + 

v    J  (t-x)  «     a 

—  CO 

after  establishing  the  interesting  fact  that  this  transform 
obeys  the  semi-group  property 

(2.3)        Pb  Pa  *  =  Pa+b  k    . 

After  the  proof  of  (2.^)  Pollard's  method  of  analysis  is  given 
by  the  following  sequence  of  equations: 


'  -  ■ 
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p  f  =  P„Pa*  .  r,,^  =  i  I     (m>  »(t)  ot 


—  CO 


(t-x)c+(y+a)' 


f(s+iy)  +  f(x-iy)  =  i   [ £_ 5  + 

T  _w  (t-x-Iy)  +  a    (t-x+iy)  +a 


"2~ 5" 


■]4>(t)  dt 


f(x+iy)  -:-  f(x-iy)  -  P  f 


-±-  F 

2iri  J 


1       1 

t-x-i(y+a)  "  t-x-Hy-a) 

,     1  1 

t-;:+i(ya)  t-x+i(y+a) 


c-x-i(y-:-a)   t-x+i(y+a) 


1   /"  (a-yl^t)  dt 
T  j4   (t-x)-+(a-y)2 

CO 

11m  [f(;:-:-iy)+f(x-iy)]  -  -i   /"'  -  f(tA  %  -■     ^(x) 
y->a- 


<t>(t)  dt 


—  CO 


(t-x)  +a 


Pollard's  proof  of  the  semi-group  property  (2.8)  does  not 
depend  on  the  \se  of  harmon:.c  functions.  The  Poisson  transform, 
however,  is  most  likely  to  arise  In  connection  with  problems  in 
potential  theory,  and  a  proof  based  on  this  theory  may  be 
desirable.  Me   proceed  to  give  a  proof  of  this  kind. 

If 


f(t)  =i   F*  Ml%*  Pa  d, 

T  £   U-t)2+a2 


then 


b  aT 


IT 


b  ,     /'a  d>(e)  dg 


-co 


(t-xr+b-  _^3  U-tr+a' 


Now  the  function 
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C3 


i       ^  If  y  l    P        <1>U) 


cV 


oux,yj     =     -   /     -5 — 9   •  -    /     -5 — *  dt 

is  a  potential  function  for  y  >  0  and 

Pb  Pa  4>  -  co(x,b)  . 
We  also  see  that 

CO 


C,(X,0+)  =  I  r    ajfcii)  dj 


—00 

CO 


=  1±m  1  r  (y+a)  ?(£,)  *  t 

y->0   Ji    (e-x)2+(y+ar 

Furthermore,  the  function 

1  f   (Y+z) J>(Z)   dg 
U  _M  (?,-x)2+(y+a)2 

is  harmonic  in  the  upper  half-plane,  and  since  it  Is  equal  to 

co(x,0+)  when  y  ->  0,  it  follows  from  the  appropriate  uniqueness 

theorem  that 


CO 


(  ,  __  1  r  to^tBj  ^1 : 

«i  (?,-x)2+(y+a)2 


We  conclude  from  this  that 

Pb  Pa  * .  c,(x,b)  .  I  f  (*■»>+(€  >  d|  .  p     4- 

b  a  ^  (£-x)2+(b+a)2    a+b 

The  Idea  of  this  proof  can  easily  be  extended  so  that  it 
can  be  applied  in  the  analysis  of  certain  integral  equations 
which  may  be  regarded  as  generalizations  of  Poisson's  integral 
equation.   Some  of  these  equations  are  considered  in  the  sequel. 

We  have  seen  that  the  semi-group  property  implies  the  inversion 
formula.  We  remark  in  passing  that  conversely,  it  is  not  difficult 
to  prove  that  the  semi-group  property  is  implied  by  the  inverseion 
formula. 
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III.   The  Poisson  transform  of  <j>  is  a  convolution  integral 
and  this  means  that  we  can  also  invert 


U)  =  *  J  n   \^  p 

(t-x)  +a 


"2~~ 2 

by  using  Fourier  transforms.  Let  (j>(t)  be  a  function  such  that 
we  can  assert  that  the  Fourier  transform  of  the  integral  is  equal 
to  the  product  of  the  transforms 


J(m)   =  f  elv;::  (|>(x)  dx 


J 

—  CO 


and 


.  iwx 


17 e  then  have 


where 


-a|w|    1  r    exv'"adx 
=  *  J         x2+  a2 

—  CO 


F(w)  =  e"a|w|  i(w) 


F(w)  =   ;rei::v;  f(x)  dx  . 


It  follows  that  the  transform  of  <J>  is 

J(\i)   =   ea'u'  P(w) 

.  [ealwl  +  e-a'v;l  -  e-a'wh  F(w) 
.  [eaw  +  e-aw  -  e-a'w|]  F(w) 

and  hence  by  inversion 

CO 

(2.9)  *(::)   -   f(x+ia)   +  f(x-ia)    -\   f    a   £{t\  Ci% 

T  (t-x)"+a- 

-co 


. 


' 


■ 
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This  proof  is  formally  simple  and  elegant j  but  the  justifica- 
tion of  the  formal  procedure  involves  a  much  more  complicated 
analysis  in  which  <j>(x)  must  apparently  be  subjected  to  severer 
restrictions  than  the  mild  ones  stipulated  in  the  previous  methods. 
For  example,  the  above  theory  at  least  requires  the  existence  of 
the  transform  ]f> (") ,   but  the  non-existence  of  jjj(w)  does  not  inter- 
fere with  the  proofs  I  and  II. 

The  formula  (2.9)  is  often  called  Snow's  formula.   C.  Snow  [2] 
derived  it  in  1924  by  using  operational  procedures  based  on 
ordinary  Fourier  transform  theory,  and  subject  to  rather  severe 
restrictions . 

The  confining  conditions  imposed  by  the  ordinary  Fourier 
transform  theory  can  of  course  be  considerably  lightened  if  we 
think  of  this  theory  in  terms  of  generalized  functions  and  their 
transforms.  The  use  of  the  more  sophisticated  theory  also  allows 
us  to  retain  formal  simplicity  and  elegance;  but  it  should  be 
remarked  that  the  theory  of  generalized  functions  introduces  an 
apparatus  which  is  much  more  advanced  than  the  methods  of  proof 
presented  here,  in  this  section.  Nevertheless,  the  results  of  an 
appeal  to  the  theory  of  generalized  functions  are  worth  noting. 
According  to  this  theory  the  inverse  of 


]5(w)  =  ealv;|  F(w) 


is 

CO  00 

«j,(x)  =^_|e-ixv'  ea|w|  feltv   f(t)  dt  dw 

-00  -00 

or 


!  . 

i 

■ 

-  -  . 


' 


(t-x)a+a2 


11 

CO  CO 

(2.10)  <|)(x)   =  I  Re  /eav'  e"1™   f  eltM  f(t)    dt  dw 

0  -00 

where  the  transforms  are  to  be  Interpreted  as  generalized  transforms 
Thus,  if  we  write 

00 

f(x)  =  pfl  d»  =  ^  r  ^-^tj  dt 

a 
and  symbolically 

a 

we  can  interpret  [P.]"  to  mean 

oo  oo 

(2.11)  1    f  _  1  Rerea«  e-ixw      r    itw  f(t)    dt  dw 

0  -co 

where  vie  assume  a  _>  0.   From  the  present  point  of  view,  the 

operator  P  f  is  defined  by 
a 

00 


e  I  Re    /P    e-aTJ  e-lxw      f  eltw  f(t)    dt   dw 
v        J  J 


0  -co 

IV.   K.  Bateman  [?]    showed  that  the  inversion  formula  for 
the  Poisson  transform  can  be  found  by  using  the  theory  of  continua- 
tion of  potential  functions.  The  following  is  a  somewhat  amplified 
account  of  his  method. 

If  a  solution  of 


(t-x)  +a' 

exists,  then  f(z)  must  be  analytic  for  -a  <   z  <  a.  The  real  part 


■ 


. 


■-'• 


1     . 
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of  f(z),  namely,  Re  f(x+iy)  is  harmonic  in  -a  <  y  <  a  and  it  reduces 
to  f(x)  when  y  =  0.   As  v;e  know,  Poisson's  equation  may  be  regarded 
as  prescribing  the  value  of  the  potential  function 

*<x,y)  .  i  /   ^Mi^       y ,  0 

at  y  =  a,  that  is  #(x,a)  =  f (x) .  How  #(x,y)  is  the  real  part  of 
a  function  P(x+iy)  which  is  analytic  in  z  =  x+iy  for  y  >  0.   In 
effect,  the  integral  equation  poses  the  problem  of  deducing 

lim  Re  P(x+ie) 

e->0 

from 

Re  P(x+ia)  =  f(x)  . 

The  function  ReP(x+iy+ia)  is  harmonic  in  the  strip  -a  <  y  <  a. 
The  difference 

Re  (F(x+iy+ia)  -  f(x+iy)} 

is  harmonic  in  -a  <  y  <  a  and  it  vanishes  for  y  =  0.  Therefore, 
by  harmonic  continuation  from  the  strip  0  <  y  <  a  into  the  strip 
-a  <  y  <   0  we  have 

-  Re  (F(x+iy+ia)-f(x+iy))  =  -   { F ( x+iy+ia ) - f (x+iy) ) 
or 

He  {P(x-iy+ia)+P(x+iy+ia)J  =  f (x+iy)+f (x-iy)  . 

If  ue   let  y  ->  a-  we  find 


■ 


- 

..  . 


.   ■■ 


• 


. 


' 


. 
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lim     Re  F(x-iy+ia)   =  ?//(x,0+)   =   <|>(x) 
y->a- 


=   lim  (f(x+iy)+f(x-iy))    -Re  F(x+2ia) 

y->a- 


Por  y  >  a 


CO 


Re  F(x+iy)   =  *(x,y)   =  \     f    [y'a]^[t]    ?| 


(t-x)~+(y-a) 
and   from  this 

CO 

Re  F(x+2Ia)   =i      f    a   f{tl  %   . 

—  CO 

Thus  once  again,   v/e  see  that 

CO 

<|>(x)   =   lim  [f(x+iy)    +  f(x-iy)}    -  i      T    a   f(tj)  d£ 
y->a-  u  j£      (t-x)2+a^ 

At  first  sight  it  may  seem  that  we  can  start  with 

Re  F(x-i-ia)  =  f(x) 
and  then  immediately  find  <j>(x)  from 

<|>(x)  =  lim  Re  F(x+is)  =  lim  f(x-ia+ie)  . 

£->0  £->0 

However,  this  is  only  effective  if  we  can  see  how  to  express 
f(x)  in  the  form 

f(x)  =  Re  F(x+iy) 

y=a 

where  F(x+iy)  is  analytic  in  z  =  x+iy  for  y  >  0.   [This,  in 

reality,  is  the  essence  of  the  problem  of  solving  the  integral 

equation;  it  is  what  the  inversion  formula  accomplishes.]   For 

example,  if 


• 


. 


■ 
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-w  x    m  1      T    a  b(t)   dt 

x  ( x )  =  s  in  x  =  —   /  — T   A — 5 

(t-x)N-a~ 


J 

—  CO 


then  the  solution  <f>(t)  is  not 


lim  sin  (x-ia+ie)  =  sin  (x-ia) 
s->0 

because  this  is  not  real  and  the  inversion  formula  shows  it  must 

be  when  f(x)  is  real.   On  the  other  hand,  if  we  recognize  that 


f(x)  =  sin  x  =  -Re  i  ea  e1::-a  =  -  ea  Re  i  ei(x+iy) 


y=a 


then 


f(x-ia)   =  -eaRe  i  eKx-ia+iy) 


=  -ea  Re  i  eix  =  ea  sin  x 


y=a 


does  In  fact  give  the  solution. 

In  connection  with  the  remarks  in  the  last  paragraph,  it  should 
be  noted  that  the  expression  of  f(x)  in  the  form 


f(x)  =  Re  F(x+iy) 


y=a 


is  immediate  If  we  use  the  theory  of  generalized  functions 
This  theory  gives 


00 


f(x)  =^Re  fe1^  J V1^  f(t)  dt  d£ 

0        -CO 

-  I  Re  |e^(^iy)  eae  J  e^t   f(.  }  flt  ^ 
0  -«> 


y=a+ 


provided  the  generalized  transforms  exist.  We  suppose  that  f(z) 
is  analytic  in  -a  <   z  <  a  and  such  that 

CO  CO 

F(x+iy)  =  I  re*(x+ly)    e*i    fe'1^   f(t)  dt  de; 


-CI 


- 


... 


-    • 

.   - 
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is   analytic   for  y  >  0.      Then 

CO  00 

<|>(x)   =   f(x-ia)   =^RejTea^    e'^x   F  eiet   f(t)    dt  d£ 

0  -co 

which  is  the  formula  (2.10).   If  f(t)  =  sin  bt,  the  generalized 
trans form  is 

C3 

F    e1^t   sin  bt  dt  =  iTr[6(£-b)  -  5(*+b)] 

—  C3 

where  5(t-t, )  is  the  generalized  function  defined  by 

00 

F  5(t-t1)  A  (t)  dT  =  A(t1)  . 


Hence 


kix)   =  f(x-ia)  =Re  i  F  ea^  e"lfjX[  5U-b)-5(£+b)  ]  de 


n  .,   ab  -ibx    ab   .  , 
=  Re  i  e   e     =  e   sin  bx 


The  above  methods  suggest  that  we  can  extend  our  considera- 
tions to  geometric  configurations  other  than  the  half-plane. 
They  also  suggest  that  problems  similar  to  the  one  discussed 
above  can  be  solved  for  other  elliptic  partial  differential 
equations  besides  the  two-dimensional  potential  equation. 


- 

( 

- 

■ 

■ 

: 

■ 
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3 .   Poisson's  Integral  Equation  for  the  Circle. 

The  equation 

p   p  '-^  c    <  v 

^0  r  +a  -2ar  cos(9-a)     -   - 

can  be  referred  to  as  Poisson's  integral  equation  for  the  circle. 

It  can  be  identified  as  the  one  which  needs  to  be  solved  in 

connection  with  the  problem  of  finding  a  potential  function 

V(p,9)  which  satisfies  p  ^oo  +  pyw  +  7qQ  =  0  for  p  <  r  and  which 

is  equal  to  f(9)  along  the  circle  o  =  a  <  r.  We  assume  that 

f(a)  is  merely  integrable.   If  iMr-,9)  =  <f>(9)  almost  everywhere 

we  have 

p   p   27T 

(3.2)        7(p,9)  =  'p_    /  — p- — p— r 

'    -q  r_+p  -2pr  cos(9-a) 

and  if  ^(a,9)  =  f(9)  then  <j>(«)  must  satisfy  (3.1). 

The  integral  in  (3.1)  can  be  written 

i---a^   f  tj>(cQ  da 

27r   ^  [a-re^^Ha-re1^] 

i9 
or,  if  we  replace  e   by  z,    as 

o   o   2tt 
r  -a    P      ft  (a)  da     

2Tr   -0  [«.r«-ao«][a.J2S] 
This  integral  defines  a  function  of  the  complex  variable 
z  =  pei6,  analytic  in  the  domain  a/r  <  |z|  <  r/a.   It  shows 
that  (3.1)  possesses  a  solution  only  if  the  prescribed  function 
f(9)  is  an  analytic  periodic  function  of  9  which  allows  continua- 
tion into  the  annular  domain  defined  above.   If  this  condition 
is  satisfied,  the  equation  (3.1)  can  presumably  be  solved  by 
using  any 


- 


■ 


: 
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using  any  one  of  the  methods  of  Section  2.   [In  order  to  solve 
(5.1)  by  the  transform  method  III,  this  method  must  of  course 
be  changed  to  a  finite  transform  method.]   Let  us  choose  method  II 
v;hich  depends  on  the  existence  of  a  semi-group  property  of  the 
transform  because  of  the  four  methods  this  one  seems  the  most 
challenging,  and  most  likely  to  involve  a  nev;  Idea. 

If  we  put  A  :  a/r  <  1  vie  can  virite  (J-.l)  In  the  form 

(j)(q)  da 


(3.3) 


2  2v 

r(9)  =iijr 


—  -   ?\  *  • 

^0   l+A  -2A  cos(a-O)      A 


With  0  <  A-,  <  1  and  0  <  A2  <  1,  the  repeated  transform 


PA  P^  <|>  is 
*1  *2 


\   A2 


The  function 


Kir2 


2ir 


2tt 

r  1 

0     i+^i-^x003^-9)    0   l+^p-2^2cos(a-p) 


2tt 

,p_Ma)dads 


J 


„(„.»,.  iigj  / 


^  1+p  -2p  cos(p-Q) 
is  a  potential  function  for  0  <  1  and 


2,    2,/ 

i 

~2w      I  77T2 


(l-^o) 


Mi 


a)    da 


•£   l+A--2A2cosfct-^) 


dp 


%%*  =  ffl(>l'Q) 


l/e  also  see  that 


2      2tt 


cod-,9)    =     -^    J 


6(g)    da 


,   l+A|-2A2cos(a-9) 


2t 


•^     /   <*>(<*)   Re   ( 


1+Ve 


i(9-a) 


1-A2e 


i(O-a) 


}    da 


- 


■ 


) 


• 


■ 
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1+T^o  e1'9"^ 


=  11m  4-  b(a)   Re  f £ if-A.^3  da 

i   p          1+ZAp  e 
=  lim  4=  /  <t>(a)  Re  ( _■  ■}  da 


;™l^iT l-z^e"1 


10 

where  z  =  o  e    .   Since 


-I    p         1+za0  e 

is  a  potential  function  for  |z|  <  1,  and  since  it  is  equal  to 
o>(l-,9)  when  jz|  ->  1  it  follows  that 

o(p,9)  =  -i  /  <><<*)  Re  { i-^)  da  . 

O  1— ZA^e 

We  can  now  see  that 

1+>   x„  ei(e^} 


,        p  j.-ta,  A-e 

P      P      $  =  co(a    ,0)    =  4-    /    <|>{a)   Re   { 1/Q-a)) 

*1   A2  -1  27r^'  l-A,7v0e119KX) 


da 
l"2e 


2,    2tt 


[1-(A1A2)    ]     p   <{)(a)    da 


J      l+(A1?v2)--2A1A2cos(a-9) 
and  hence   that  P,(b  obeys   the   semi-group  property 

A 

Al  A2     Al*2 
This  property  can  be  used  to  invert  (J.l)    in  the  following 
way.   The  equation  (3".D  can  be  written 

1    2p  ia  .   -i9 


Then 


• 


■ 


. 


i 


■ 
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f(0+i  In  o)    +  f(©-i  In  p)   -  P  f  =   f(0+i  In  p)+f(9-i  In  p  )-PpA4> 


f(9+i   In  p)    +  f(9-i   In  p)    -  P   f 


2tt 

n 


'   2itJ 
0 


<-.  e" 


ie 


+ 


Ape 


-19 


peia-Aei0 

eia 

+     ia   .       19        + 
e     -Ape 

ia 
e 

e      -Aoe 


e~ 

ia    ,       -19 
-Aoe 

.       -19 
A   e 

pe 

Ace"19 

.        .1            .    .          .      -       .    _ 

e        -Ape 


Ma)    da 


from  which 


2?r 


f(0+i   In  o)    +  f(9-i   In  o)    -    (1;£    }    f    «-£ 


(a)    da 


_    (p    -A    ) 
2tt 


^7T 


1+p   -2o   cos(a-9) 
<}>(a)    da 


— 5 — 5 
JQ     p    +A   -2pA   cos(a-9) 


Finally,    if  we  let  p    ->   A+  we   find 


(3.4) 


<j>(9)   =   lim     (f(9+i  In  p)+f(9-i  In  p)) 
p->A+ 

o        2tt 
_    U-*    )       f  f(a)    da 

4,   1+A2-2A   cos(a-9) 


for  almost  all   values   of  9.      This   is   the   inversion  formula   for 


(3.3). 


' 


. 


- 


.    • 


• 


■ 


?)' 


■ 
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■ 

- 
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4 .   Some  Generalizations  of  Poisson's  Integral  Ecuation. 

Poisson's  Integral  equation  suggests  a  more  general  problem. 
Let  D  be  a  simply  connected  domain  whose  piecewise  smooth  boundary 
C,  is  given  by 

z  =  z,  (t)  a  _<  t  <  b 

and  let  Cp  ,      given  by 

z  =  Zp(-r)  c  <   t  <  d  , 

be  a  simple  piecewise  smooth  closed  path  in  the  interior  of  D. 

Let  V(x,y)  be  harmonic  in  D.   Suppose  that  at  the  boundary  C, 

the  limit  of  #(x,y)  is  (j)(t)  and  that  if  (x,y)    is  a  point  in  D 

then  Tp{x,y)    is  given  by 

b 
f(x,y)   =JG(x,y;t)  <j)(t)  dt  . 
a 
The  problem  we  have  in  mind  is  this.   Solve  the  integral  equation 

b 
(4.1)   f(T)  ■  Re{H[z2(T))  =  J   G[x2(t)  ,y2(x) ;  t]  <|>(t)  dt    c  <  t  <  d 

a 
where  z„(t)  =  ;cq(t)  +  iy„(r)  is  on  Cp.   We  assume  that  H(z)  is 

analytic  for  z  on  Cp  and  for  z  in  the  domain  bounded  by  C?  and  C,  . 

We  also  assume  of  course  that  the  prescribed  function 

f(f)  b  Re(H[z0(f)]]  possesses  all  of  the  properties  implied  by 

the  integral  representation  on  the  right  hand  side  of  the  equation. 

If  we  set 

f(x,y)   m   Re(P(z)) 

where  z  is  an  analytic  function  of  z  =  x+iy   for  z  in  D,  then  the 
solution  of  the  integral  equation  (4.1)  requires  the  deduction  of 


. 


" 


-     ' 
.  .  ... 


■ 


' 


. 


t 


- 


' 


: 
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11m   Re(F(z))  =  <j)(t) 
z-> 

from 


z->z, ( t) 


Re(P[z2(T)]]  =  f(-r)  =  Re{F:[z2(x)]}  . 

This  can  be  accomplished  with  the  aid  of  conformal  mapping. 

Let  £  =  m(z)  be  the  function  which  maps  the  domain  enclosed 
by  Co  into  the  upper  half  U  of  a  £ -plane.   Under  tl is  mapping  the 
image  of  C,  will  be  C  ,  in  general  a  curve  not  in  U,  but  in  the 
Riemann  surface  formed  by  properly  connecting  superimposed 
£ -plane  sheets  in  accordance  with  £  =  m(z).   However,  v;e  assume 
here  for  simplicity  that  C  is  in  the  lower  half-plane  on  the 
same  sheet  as  U.   The  image  of  D  is  then  the  domain  D*  which 
contains  the  point  at  infinity  in  the  £ -plane  and  which  has  C* 
for  its  boundary. 


Let  the  inverse  of 


be 


C,   -  m(z) 


2  =  M(0 


With  this  we  have 

Re{F[MU))    =  Re{K[M(e)]]    . 
Now  the  difference 

Re[F[M(01    -  H[M(?)]] 

is  a  harmonic  function  in  the  part  of  the  exterior  of  C*  which 
lies  in  the  lower  half  plane;  and  this  function  vanishes  for 


- 
■ 


■ 


■ 


• 


■ 


' 


■ 
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y 


\ 


—  _Hi 


z   = 


z2(t) 


^ 


z-plane 


z  =   z1(t) 


-> 

x 


Ug.    1 


Tl 


£ -plane 


-&■ 


/' 


\ 


.    * 


Fis.   2 


■ 


<■ 
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Im  £  =  0.   Kence  if  t,      is  a  point  on  C  ,  harmonic  continuation 
across  Im  t,   =  0  gives 

11m     Re{P[M(0]    -  H[M(C)]}   =   -   11m     Re(F[M(C )  ]-H[MU)  ] ) 
or 

lim  4lRe[F[M(0])   =   Hm  ,,Re{H[M(£)  ]+H[M(C)  ] )    -  Re(P[M(^*)]} 

c->c  c->e 

For  Im  £  >  0,  Re{F[M(£!)]}  is  harmonic  and  its  value  at  Ira  ^  =  0 
is  Re(H[M((;)]].   Therefore 

-i       i  - 1 

A  return  to  the  z-plane  gives 

(4.2)        «|>(t)  =  lim    Re(F(z)) 

z->z-,  (t) 


=   lim   Re{H(z)+K(M[m(z)])} 
z->z1(t) 

,     rd  Re{K[z„(T)]}d  m[z„(x)] 
-  -  Im  i g  — 

71    Jc    m[z2(T)]-m[z  (t)] 


or 


(4.3)        <|>U)  =  lim    Re[I!(z)+H(M[m(z)])} 

z->z, (t) 


-1  Im  J  Re(R[z2(T)]}  ^|  ln{m[  z2(t)  }-m[  z^t)])  dr 
c 
for  the  solution  of  (4.1). 

H.  Bateman  [  3]  considered  the  special  case  In  which  C2  is 

given  by  z  =  z2(t)  =  t  +  in;  C,  is  given  by  z  =  z,(t)  with 

Im  z, (t)  <  h;  and  D  contains  the  point  at  infinity  as  shown 

in  Fig.  3- 


. 


' 


. 


' 


' 


• 
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t 


r 

h 


D 

/" 

"s. 

• 

1 

e* 

Fig.  5 


Here  the  mapping  function  is 


£  =  m(z)  =  z  -  ih 
z  =  M(C)  =  C  +  ih 

and  (4.2)  shows  that  if 

b 
#(x,y)  =  j/  G(x,y;t)  <|>(t)  dt 

a 
is  a  regular  potential  function  in  D  with  boundary*-  value  <}>(t) 

on  C-j,  then  the  solution  of 

b 

f(T)  =  f[z0(T)-ih]  =  f     G(T,h;t)  <t>(t)  dt 


is 


t(t)  =  lim     (f(z-lh)+f(5+lh))  -  i  im  f  ^jll^V)-;  .i 

This  agrees  with  Bateman's  result. 


. 


• 
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The  above  results  can  be  extended  and  specialized  in  various 
ways  but  we  will  consider  only  one  more  simple  specialization 
which  may  be  or  interest.  Suppose  that  a  part  r  of  the  curve  C0 
is  brought  Into  coincidence  with  a  part  of  C,  ,  and  that  the  remain- 
ing part  of  CD  becomes  a  line  segment  L.   For  this  case  we  can 
propose  finding  the  value  of  #(x,y)  on  C,-T,  where  ip   is  a  potential 
function  in  D  with  its  values  prescribed  on  r  and  L.   The  method 
of  this  section  can  still  be  applied  and  the  first  step  is  to  map 
the  domain  bounded  by  r+L  into  an  upper  half  plane.  By  way  of 
illustration  let  us  find  ^(x,0)  =  <j>(x)  for  x  >  0  if  #(x,y)  is 
known  to  satisfy 


>2 


>2„ 


-  +  °  y  -  o 
2    T?-° 


y  >   0 


tf(x,0)    =   0 


<   0 


i>{r  cos  c. ,    r   sin  a)   =   f(r)        r>0    ;  a    <  -S 


lT 


Fig.    4 


crr 


L 


a 

\ 


' 


• 
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The  function  Vy(x,y)  is  given  by 


"  J   (5-x)2+y2 


and  If  zp(f)  =  re-  ,  the  problem  is  to  solve 


f(T)  =  f[e-toz0(T)]  =  1      f   ±E  sing)  <j>(«)  df, 


7T        „2,  2  r 

-f,   £  +r  -2%x   cos  a 

00 

i    r  kg)  dg 

0   '' 


for  <{>(!;)•  The  function  f(w)  must  be  analytic  for  w  =  t  and  as  we 
can  see  from  the  integrals,  f(r)  must  allow  analytic  continuation 
to  f(w)  where  vj  is  complex  and  -a  <  arg  w  <  a.   The  function 


o 


C  =  [z  e"10]7'"-0  =  m(z) 

Tr-a 
2  =  eia  C  1T     =   M(C) 

maps  the  section  bounded  by  r+L  into  the  upper  half  of  a  £ -plane 
sheet  and  under  this  mapping  the  image  of  z  =  p,  o  >  0,    is  the 
ray  £  =  p  exp{  -Icir/ir-a)  in  the  lower  half  of  the  same  sheet.  With 

/  %      ia 
z2(t)  =  t  e 

z-^(t)  =  t 

H(z)   =  f(z  e_la) 

■ 

the  formula  (4.2),  that  is, 

4>(t)  =  lim    Re[K(z)+K(M[m(z)])) 
z->z, (t)  , 

1       -     rd  Re(r[zo(T)]]  d  m[zQM] 
-  -  Im 


71"     J 


m[z2(r)]  -  m[z1(t)] 


. 


■ 


■ 


. 


• 
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Is  applicable;    and  we   find 


<|>(t)   =   11m  [f(t   e"i9)    +  f(t   eiG)] 


9->a  „  a 

Im 


-i  n  *.•     \       7T-Q     3 

1       -        r  f(T)     T  dT 


7c-a 


0       ir-a        r .       ia  n  -rr-a 


The  reader  will  perhaps  find  it  interesting  to  supply  the 
modifications  which  are  necessary  when  -5  <  a  <  it  . 


5.   The  Boundary  Values  of  a  Biharmonic  Function  which  is  Prescribed 
Along  Two  Interior  Paths. 

Problems,  analogous  to  those  studied  in  the  previous  sections, 

can  also  be  proposed  for  the  two-dimensional  biharmonic  equation 

which  is  a  dominant  equation  in  the  theory  of  elasticity.  We 

proceed  to  discuss  one  of  these  problems  and  thereby  suggest  some 

other  applications  of  the  inversion  formulas  found  above.   The 

problem  we  choose  is  this.   Show  how  to  find,  in  closed  form, 

formulas  for  :;[r,Q)   and  #  (r,Q)  if  it  is  given  that  :!/(?, 0)  satisfies 

in    >    i-n2    i-n    5s,    i-\2,         0<o<r 
fl  o    d  ,  1  o  1  r  1  o    o*  .  1  3  -!/,   r,  —  ■ 

7//(a,9)  =  f(0)  0  <  a  <  b 

-,'/(b,9)  =  g(6)  a  <  b  <  r  . 

If  i'/[o,9)    is  a  biharmonic  function  in  0  _<  p  <   r  it  is  well 
known  that 


■ 


'     - 

- 


■ 


.   ' 
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,  2  2.2  ^  [r-o  cos(a-Q)]  ^(r,a) 

*(p,©>  -  (%Zg  }    /  — r 5-g  dec 

^     ^0  [r  -2rp  cos(a-9)+o  ] 

(r2-o2)2  f"    V*V*>  ** 

^7rr    ^  r2-2-'ro  cos(a-9)+r.2 

Hence  the  problem  just  proposed  is  equivalent  to  solving  the 


following  pair  of  integral  equations 

La2j2  ^[r-a   cos  (a -9)]    tf 
2irr     ./  77T7~T~7TT77 


^   [r   -2ra   cos(a-9)+a"~] 


2?r 


(r2-a2)2      f        Vr'a)    dci 
•i         zprr; /   — 73 ' 3 

^  r~-2ra   cos(ci-9)+a~ 
,2  .2,2     2J  [r-b   cos (a -9)]    ^(r,a)    da 

g(9)  =  (rc; }      /  —g ^—  + 

^     [r  -2rb   cos(a-9)+b~] 

.,.    (r2-b2)2     f  ^,a)    *»  > 

^7rr  ■£     r2-2Vo   cos(a-9)+b2 

These   equations   could  be   solved  by  using   finite   transform 

techniques;    but   the   course   to   the   solution  is   easier  and  more 

satisfactory  if  we  use  the  fact  that 

#(x,y)   -  Re  [z  F^z)    +  G1(2s)3 

is  a  necessary  and  sufficient  condition  for  ^(x,y)  to  be  a 
biharmonic  function  in  a  domain  D  provided  that  F, (2)  and  G(z) 
are  analytic  for  z  =  x  +  iy  in  D.  The  last  equation  can  be 
rewritten  as 

V(::,;:)  ■  (r  -zz)  Re  F(z)  +  Re  G(z)  , 

a  form  which  is  especially  appropriate  for  the  analysis  of 
functions  which  are  biharmonic  in  a  disc  with  center  at  the 


. 


- 


• 


. 


29 

origin  of  the  coordinate  system.   For  our  purposes  it  is  most 
convenient  if  we  take 


tip ,9)   =  (r2-o2)  ne-^L-  f   -^jjjif}  n(^  In  |)  d^ 


(5.D 


v.'here 


r 
c: 


+  Re  ik  f  xxt^  v  (i ln  I'  ^ 


,,      ia  i© 

£  =  r  e    ;  z=pe 

and  each  of     \i[-r  In  •£)    =  p.  (a)    ;    v  (-r  In -^)   =   v(a)    is  real.      From 

(5.1)  we  find 

(5.2)  *p(p#0)   -   -2?  Re  -i   •  £    ^jj+|>     ^(1  m  |)   d£ 

irp  d©         J"  1,-z 

<r~  Re  ±r  r, . 

iro    o©         j  l,   -   z 

C 
If  we  let  p   ->  r  in   (5.1)   and   (5.2)   we  obtain 

^(r,9)   =  v(Q) 

f0{v,9)   =   -  2r  n(©)   -  i  Re  <£   v(%lgdg 

c 

which  show  how  ^(r,©)  and  Vy  (**>©)  depend  on  u.(©)  and  v(©)  . 
The  problem  has  now  been  reduced  to  solving  the  system 


' 


■ 
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f(e)   .    (r2.a2j   Re^      f     (^e-)^)dC+Re^    fJL±«      vWflC 


g(0)   =    (r2-b2)   Re     1      ri^e19)a(a)    dC   +  Re  ^    f     <  ^^    v(a)dC 


777—10; —  +Re  WlJ 

C  ^U-be      )  c 

where   £  =   r  e      .      These   equations   reduce  to 

, _2     2*      ^     [(r2-a2)  n(o)+v(o)] 

f (0)   =    (r  "a    }      /     g 5-  da 

dir  ^         r  -24a  cos(a-0)-:-a 

,2     2,      2rT  [(r2-b2)   n(a)+v(a)] 
g(e)  =  (r  g  }     /  ?— - —2-   <*>   . 

^ir         ^  r   -24b   cos(a-e)+a 

They  can  be  solved  by  using  the  inversion  formula  given  in 
Section  *.  Me   find 

lira     [f(0+  i  ln-g)  +  f(©-  i  ln£)] 

o->r- 


(r2-a2)u.(9)+v(0) 


(r2-b2):-i.(9)+v(0)    - 


P      p      2tt 
(r^-a-)    T  f(a-)    da 

-    nT- '         — 2 ~~~^>     > 

d"         JQ     r   -  2ra  cos(a-Q)-:-a£ 

11m     [g(©)  +  i  ln|)  +  g(©-i  ln-g)] 

o->r- 


P      o        27r 
(r   -b^)       f  g(q)    da  . 

27r         J      r2-2rb   cos(a-9)-:-b2 


and  from  these 


• 


- 


(b2-a2)v(9) 


lim 

o->v~ 


lvd-ef)    [g(9+i  ln|)   +g(©-i  ln|)]     J 

c 

I   -    (r2-b2)[f(0+I  ln|)    +  f(0-i  ln|)]  \ 

(r2-a2)(r2-b2)       P  f(a)    da 

__ — _ _ 


zir 


i  r2-24a   cos(a-9)+a' 
2ir 


(r2-a2)(r2-b2)      P  g(a)   da 


2v 


0 


_ ^ 

r   -24b   cos(a-e)+b_ 


(b2-a2)n(9) 


i'     f  (0  -i-  i  In  £)   +  f  (9  -  i  In  £)    7 

j  a  a     i 

lim        x 

\   -g(9  +  i  ln|)    -  g(9   -   i  ln-g) 

V  / 

o      P        2tt 
,    (r   -b~)       P  E(a)    da 

+ ST"    .        .2   .._,..    ^..2 


p      p 


o 
2tt 


r   -2rb    cos(a-9)4-b 


(r~-a~)       /"    f(a)    da 


J     ~5~ 


r  -2ra  cos(a-9)+a' 


■ 


■ 


- 
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6.   The  Inversion  Formulas  in  Forms  Suitable  for  Applications. 

'..re  have  used  several  different  methods  to  show 

(a)  if  <|>(t)  is  integralle  over  each  finite  interval 

(b)  if  each  of  '    ${tl   dt  j   /  _£i_y_di_  exists  with  c  >  0 

-co       t  co      t 

then  the  solution  of 

C3 

(6.1)  f(z)  -  |  f  iiispi  ""  *  *  *  " 

u  J   (t-x)N-a-  a  >  ° 

—  CD     v 

is 

CO 

(6.2)  M-)  =  lim  [f(x+iy)  +f(x-iy)]  - -|  f  a  r(U  d| 

y->a-  F  _4  (t-::)2+a2 

for  almost  all  values  of  x  provided  f(z)  is  analytic  for 

-a  <  Im  z  <  a.  The  "almost  all"  stipulation  can  be  removed  by 

contracting  the  class  of  functions  to  which  (f>  may  belong*  Let 

us  assume  hereafter  that  instead  of  the  condition  (a)  the  function 

(|>(t)  is  required  to  satisfy  a  uniform  KSlder  condition  for  all  tj 

and  that  in  addition  the  prescribed  function  is  such  that 

lim  f(x+iy)  =  f(x  +  ia)    and     lim  f(x-ly)  =  f(x-  ia)  . 
y->a-  y->a- 

Under  these  conditions,  and  with  condition  (b),  the  inversion 

formula 

00 

(6.3)  <J>(x)  =  f(x+ia)  +  f(x-ia)  -^  f  a  X-(H  % 

—  00      ' 

holds  for  all  values  of  x. 

V/ith  respect  to  the  harmonic  function 

C3 

(6.4)        fix,7)  =\    f  y  4><^  % 

_„   (t-x)2+y2 


- 


. 


• 


• 


' 


. 


! 


; 


which  is  harmonic  in  the  half-plane  y  >  Oj  and  uniquely  determined 
in  y  >  a  by  -,!/(x,a)  =  f(x);  that  is,  by 

(6.5)       *(x,y)  -±   r3l7za)/(t)  d| 

_i  (t-x)^-(y-a)2 

the  inversion  formula  (6. 5)  is  a  continuation  formula  which  allov;s 
us  to  continue  ^(x,y)  from  the  half-plane  -co  <  x  <  ^  ;  y  _>  a  into 
the  strip  -«  <  x  <  »;  0  <  y  <   a  . 

The  inversion  formula,  or  continuation  formula  (6. 3)    is  quite 
satisfactory  if  f(x)  is  analytically  prescribed  for  -co  <  x  <  ™-, 
but  if  f(x)  is  a  tabular  function  derived  from  observations,  the 
formula  cannot  be  used  without  modification  because  f(x+ia),  for 
example,  is  not  determinable  from  a  finite  number  of  real  values 
of  f (x) .   In  addition  to  this  there  is  also  an  inherent  difficulty. 
This  second  difficulty  manifests  itself  when  we  realize  that  a 
small  change  in  f(x)  may  correspond  to  a  large  change  in  <{>.  For 
example,  if  <|>(t)  =  cos  nt  then 

CO 

1  r    a.cos  nt  dt  =  -an  CQS  ^ 
17    J      (t-xP  +  a2 

—  CO     x 

and  hence  although  the  amplitude  of  variation  of  cos  nt  does  not 
change  with  n,  f(x)  =  e    cos  nx  may  be  made  as  small  as  we  please 
by  taking  n  sufficiently  large.  We  therefore  must  admit  that  a 
small  error  in  the  measurement  of  f(x)  does  not  necessarily  imply 
a  small  error  in  the  computation  of  <t>.  This  kind  of  instability, 
for  integral  equations  of  the  first  kind  with  difference  kernels, 
was  investigated  by  Kreisel  [  4  ]  who  met  the  difficulty  by 
replacing  the  unknown  function  with  a  smoother  function. 
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Our  purpose  here  is  to  show  how  to  derive  some  modifications 
of  (6.;;)  under  the  assumption  that  errors  in  the  data  on  f(x)  can 
be  disregarded,  or  that  rapidly  changing  solutions  are  ruled  out. 
(If  a  smoothing  process  is  necessary  it  can  be  applied  after  a 
modification  of  (6.5)  has  been  obtained.)  We  do  not  propose  here 
to  enter  an  extensive  numerical  analysis  or  error  analysis  of  any 
pertinent  formula;  our  aim  is  rather  to  suggest  elementary  methods 
which  may  lead  to  practically  useful  formulas  which  can  be  substi- 
tuted for  (6.3;).  Ue  also  assume  that  a,  the  distance  from  the 
datum  line  to  the  boundary,  is  known.   If  this  is  not  the  case, 
that  is  if  f{x.,y)    is  actually  harmonic  in  a  domain  which  contains 
the  half-plane  y  >  0  ;  say  the  half-plane  y  >  -h  (h  >  0  )  ,  then 
(6.3)  and  its  modifications  simply  give  downward  continuation 
formulas  which  enable  us  to  establish  the  new  datum  line  y  =  0  . 

There  are  at  least  three  ideas  which  can  be  used  as  guides 
to  find  approximations  for  (6. J): 

(a)  find  a  suitable  approximation  for  f (x+ia)+f (x-ia)=2  Re  f(x+ia); 

(b)  find  an  efficient  formula  which  allows  extrapolation  from  the 
known  values  #(x,y),  y  >  a  ,  to  the  unknown  <j)(x)  =  <'/(x,0+)  ; 

(c)  use  or  design  iteration  procedures  which  can  be  effectively 
applied  to  Poisson  integral  equations.  These  ideas  are  not 
independent.  Each  can  in  fact  be  expressed  in  terms  of  manipula- 
tions performed  on  the  exponential  eau  which  appears  in  the 
formula  (2.10);  however,  it  seems  that  it  may  be  fruitful  to 
categorize  the  ideas  as  we  have.  Although  we  are  primarily 
interested  in  the  third  idea,  a  brief  discussion  of  the  other 
ideas  is  relevant. 


■     ' 

i 


In  the  first  place,  the  sum  f (x+ia)+f (x-ia)  =  2  Re  f(x+ia)  in 
(6.5)  can  be  replaced  by  its  Taylor  expansion  so  that  we  have 

♦  ,,,  .  2  Re  EI  -^  fn(x)  -  2  f  *  **l  % 

-CO  x  ' 

,6.6)  «,)  -  s  2Z  -^0^  f^.)  -  2  f  iittpl  . 

—  CO 

This  suggests  that  it  may  be  feasible  to  introduce  difference 
quotients.   If  we  retain  the  first  two  terms  of  the  Taylor 
expansion  in  (6.6),    that  is,  take 

CO 

(6.7)        *(x)  =  2  f(x)  -  a2  f"(x)  -  i  f  a  f(U  d| 

and  replace  f  (x)  by  the  central  difference  quotient 

r"tv\        f(::+a)-  2f(x)  +  f(x-a) 

I  (A)     «v  ^ 

a" 
we  have 

CO 

<j>(x)  =  4  f(x)  -  f(x+a)  -  f(x-a)  -  1      F    a  f(^  d| 

which  may  be  useful  when  a  is  small.   In  this  formula  we  suppose 
that  the  integral  is  to  be  found  by  a  process  of  numerical  integra- 
tion. The  last  formula  can  be  written 

(6.3)       <Kx)  =  4  f(x)  -  f(x+a)  -  f(x-a)  -  tff(x,2a) 

which  is  the  simplest  formula  for  the  downward  continuation  of 
•cV(x,y)  when  this  function  is  known  along  y  =  a  and  y  =  2a.   It 
represents,  of  course,  no  more  than  the  substitution  of  the 
simplest  partial  difference  equation  for  Laplace's  equation. 


• 
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It  may  be  possible  to  approximate  t!  e  tabulated  f (x)  with  a 
function  A(x),  say  a  rational  function  approximation,  such  that 
A(z)  is  analytic  in  a  <  Im  z  <  a.   If  this  is  the  case  we  can 
replace  f[z-i(a-s)]  in  (2.6)  by  A[z-i(a-e)]  and  use 

CO 

(6.9)  4>U)  =     "a     lm|     f  ALf^lgfLil  <^ 

£„>0    w    J  z-x-ie 

—  00 

to  compute  <|>(;:) . 

The  derivatives  of  f(x)  in  (6.6)  can  be  replaced  by 

derivatives  of  -!/(x,y)    with  respect  to  y.   We  have  f(x)  =  ^(x,a) 

and  since  ^(x,y)  is  harmonic  for  y  >  0  it  follows  that 

f"(x)  =  z'j      (x,a)  =  -  ill      (x.a)  .   Using  these  facts  we  see  t'rat 
■  xx         yy 


Om 

{    ,,111  -2m/  »  _  S   ft(x,y) 


y=a 


For  example,  (6.7)  can  be  written 


(6.10)  <|>(X)  £  2f(x)  +  a2  *_(x,a)  -if  a  f(PI  d|  . 

—  CO    *  ■ 

p 
If  a  #  (x,a)  is  replaced  by 

a2#  (x,a)  =  <>(x,3a)  -  2#(x,2a)  +  ^(x,a) 

then  (6.10)  becomes 

(6.11)  4>(x)  =  ^(x,3a)  -  5^(x,2a)  -  *^(x,a)  . 

If  A  is  the  difference  operator  defined  by  Aa>(x,y)  =  ■:  (x,y+a)-ci:(x,y) 
the  right  hand  side  of  (6.11)  can  be  written  as 

2 
?Hx,a)  -  A#(x,a)  +  A  ^(x,a) 
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and  this  is  the  same  as  P2(C)  where  P2(y)  is  the  second  degree 
polynomial,  Fg(y)  =  aQ   +  a,y  +  a0y2,  which  satisfies 

P2(Ka)  =  ^(x,Ka)       k  =   1,2,5. 

The  formula 

<})(x)  =  #(x,a)  -  A^(x,a)  +  A  V^(x,a) 

is  then  an  extrapolation  formula  from  known  values  y(x,y),  y  >  a, 
to  $>{x)    .   If  we  use  a  polynomial  of  the  n-th  degree,  that  is, 

F^(y)  =  -■   ct„y'  such  that 
n      t^s 

Pn(fCa)  =  ^(x,Ka)       K  =  1,2,  ...,n+l 
we  find 
(6.12)  <|>(x)  =  21  (-1)K  AK  *(x,a)  . 

If  we  admit  the  use  of  derivative  values  of  ?//(x,y)  for  y  >  a 
it  is  clear  that  a  large  variety  of  other  types  of  extrapolation 
formulae  can  be  devised. 

Let  us  turn  now  to  an  iteration  procedure  based  on  an 
analysis  of  Poisson*s  equation 

CO 

(6.i3)  f(x)  =\  r a  ^i  % 

TJ      (t-x)N-a2 

—  CO     x 

when  it  is  written 

Pa  ♦  -  f  . 

Symbolically, 


• 


?8 

<Mx)  =  JT"  U-PJK  f  • 

^3  a 

Consider 

(6.14)  <L  (x)   -  T1-    (l-Pa)K   f   • 

Since 

(1-P  Jf  =  #(x,a)    -  i^(x,2a)   =   -  A#(x,a) 

the  equation  (5.14)  is 

n     ft=6 
which  is  the  same  as  (6.12).  We  expect  that  if  <J>(x)  is  subject  to 

certain  restrictions  then  lim  <tn(x)  =  <|>(x)  .   If  <l>n(x)  is  known 

n->co 

then  <j>    .-i(x)    is   given  by 

n  ".1-1 

=  f  +  >-  (1-P  )  +1f 
=  f  +  a-? j  y~  ^^JKf 

a    fej  a 

or 

.     Z3     a   <t>    (t)    dt 

(6.15)  *n+i<x>  -  rM  -  *n'x'  -ij     j    ".-■>.    2  I 

—CO  *  A  ' 

a  recurrence  formula  for  successive  ^'s  •  The  last  equation 

also  shov;s  that  h     satisfies  the  Integral  equation 

n 

-,  r°°  a  (J)  (t)  dt 
^    (t-x)  +  a 

—  00 

=  f(x)  -  (l-Pa)n+1f  . 
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If  fy{t)    satisfies  (6.1;)  then  (6.l6)  is  the  same  as 


pa<+n-*>  =  -  d-pa)n+1  pa  * 


which  implies 


(6.17)  <|>n(x)  -  <|>(x)  =  -  (1-Pa)n+1  <i>  . 

If  we  use  the  property  PaPv,  =  P„,y.>  equation  (6.17)  becomes 


a  b         a+b- 

, i-(n+l)P     + 
Kn  T  'a  ni  2a 


(|>n(x)   -  <|>(x)   =   -   (l-(n+l)Pa   ++~^  PQP    ...}   ♦ 


iim  r  ^  j,  ^  -#?n?s).(<a+s!2  --J  *(t)  dt 


e  ->0  ^M   ( t-x )   +e         (t  -x)  +(a+e  ) 

00  CO 

=  -  1  lim     /"'[    re-e7v(l-e"a>v)n+1   cos  A(t-x)dA]<t>(t)dt. 

—CO  —CO 

It  can  be  shown  that  lim  <j>r(x)  =  <j)(x)  if  <|>(x)  satisfies  one  of 

n->oo 
several  different  conditions;  but  for  simplicity  we  are  now  going 

to  assume  that  <J)(x)  can  be  expressed  as 

M:0  =  P5  ^ 

where  (L  is  absolutely  inte-rable.  This  implies  that  <|>(z)  is 
analytic  in  -8  <  Im  z  <  5.  For  this  case  (6.17)  becomes 

4>n(x)  -  +  (x)  =  -  (1-Pa)n+1  Pg  4>! 
<b  (x)  -  <b(x)  -  --   f  [  f  e"5?v(l-e"aA)n+1cos  M  t-x)d*]<t>,  ( t)  dt. 


Now 


. 


• 


.. 
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oo  N  a 

yV5*(l-e-aX)n+1   cos  Mt-x)   clA  <jre-5X(l-e-aX)n+1dX+  [ 


0  a-i  kit    N 

<   (1.e-aH,n+l    d-e-°H)   . 


e  dA 


-5M 


where  N  is  some  fixed  number  independent  of  n;  and  0  <  A,  <  N  • 
The  number  N  can  be  chosen  so  large  that  e"0iV/5  is  as  small  as  we 
please,  and  then  by  choosing  n  sufficiently  large  (1-e    ) 


can  be  made  as  small  as  we  please  since  0  <  1  -  e 


-aA. 


<  1  .   It 


follows  from  this  that  <{>  (x)  converges  uniformly  to  c(>(x)  as  n  ->  «>  . 
From  equation  (6.17)  we  see  that  <j>  (x)  can  be  written 

<.„<*)  -   (n+l)Pa  -ISgiB  P8a  ...  (-«»  P(w1)bC  * 

or  since  P„<j)  =  f 


00 

4  (x)  =  (n+l)f  +i  n-i^lS ^r-g...   (-Hnn%3f(t) 


dt 


(t-x)~+a    (t-;:)~+(na)' 


(6.18)   <j>  (x)  =(n+l)  f(x) 


n 


+  i/f(t) 


—  CO 


( n-i-1 ) n  -a£    (n+l)n(n-l)„-2ar 
-  ^-j —  e    +  ^ e 

Cm  *   • 


(-l)n  e"na* 


cos£(t-x)d£dt. 


The  inner  integral  of  the  last  equation  is  equal  to 


00 

I  =  Re  f 


(n+l)n  -2a£  ,  (n+l)n(n-l)   -?af 
— ?n —  e     -i 


-=rr 


..   (-l)n  e-(n+1)a^ 


}[a+i(t-x)]«  6. 


and  two  integrations  by  parts  gives 


- 


■ 


. 
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[a+i(t-x)]i(t-x)  JQ 
Now 

CO  1 

J  a  J 

0  0 

_  i  r(i-ip)  r(n) 

a   r(n+l-ip) 

=  1  T(n) 

a   (n-ip) (n-l-ip) . . .  (1-ip) 

where  p  =  (t-x)/a.   Kence 

T  •-  Rp  f-    x     +  ( i ) n  ( n+1 )  1       ., 

1  -  Re  l   a+i(t-x)  +  a(,3-i)P(i3+i)(p+2i)...(P+ni)J  * 

The  substitution  of  this  in  (6.10)  gives 

CO 

(6.19)      <L(x)  =  (n+l)f(x)  -|  f    a  f(H  d| 


-co 

M 


(t-x)^+a' 


(n+1)  I  p-mn  r        f(x+ag)  dj3 
+  _ ir ReU)     (p-i)p(p+i)i+2i)...(3+nij 


If  we  use  the  fact  that 

CO 


(n+1)  J   PoM\n   f  d6      /„,,!,, 

—  Re(i)  J     (p-i)p(p+il...(p+ni)  =  "(n+D+2 

—  CD 

then  (6.19)  can  be  written  in  the  numerically  more  useful  form 

CO 

(6.20)      <|>    (x)   =    2   f(x)    -    /      a   r{tl  dS 


(t-x)~+a 


eo 


+  (n+1)  I  Da/^n    f 


+1)1  B,mn    /  [f(x+a;3)-f(x)]    dg 

* ReU)     J     7p-i)p(a+i)(p+2i)'-.(3+ni) 


This  is  an  approximation  which  may  be  useful  for  numerical  purposes 
if 


• 


• 


* 


• 
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(l)n 
Re  (p-i)(3(P+i)(p+2i)...(p+ni) 

is  tabulated  for  some  values  of  n. 

It  is  interesting  to  notice  that  the  approximation  (6.14), 
namely, 

n_ 

can  also  be  obtained  by  starting  with  the  solution  formula 

CO  CO 

(6.21)  <|>(x)   =  -  Re    f  eaw    /   eiv;(t-x)    f(t)   dt  dv; 


n  fen         a 


TT  u 


derived  in  Section  2.   The  exponential  eaw  can  be  expressed  as 

aw      1       ^-2-  /,   -awN'<  ,   awn   -aw»n+l 

e   =  — —  =  > U-e   )   +  e   U-e   ) 

l_(l_e"av;)   feO 

and  if  we  take 

CO  00 

(6.22)   J>n(x)  =  ±Re  f    <T"  (l-e"aV/)'^  j'  elv,{t~x)    f(t)  dt  dw 

0  -oo 

we  can  show  that  (6.22)  is  the  same  as  (6.14)  if  we  grant  that  the 
order  of  integration  can  be  changed. 

The  remarks  of  the  last  paragraph  indicate  that  various 
iteration  schemes,  or  numerical  approximation  formulas,  can  be 
developed  after  eaw  in  (6.21)  is  replaced  by  some  suitable  approxi- 
mation.  For  example  it  is  interesting  to  examine  the  consequences 
of  talcing 

Ja-e)w    (y-e)w  -(y-a)w  _  r^-  [  (y-s)w1R  -(y-a)w 

^=-  [-(y-e)]n  ^n  e-(y-a)w 

"TH Tin  e 


n=0    Ui         by 


,.„.... 

- 


.         .     •  •  .  - 


•    •■  '   - 


■ 


. 


• 


i 


. 


• 


This  suggests 

00  CO 

4>(x)  =  llm  £-   H^»n4-He  /;-<y-a)Wreiv;(t-x)f (t)  dt  dw 

~->0  ^0     n'     3yn  T    ^       _J 

or 

(6.23)  +(,)  =  11.  —UX^nHJ^l   f   (y-a)  f(t)  dt 

where  it  is  to  be  understood  that  y  >  a  .   It  is  amusing  to 
notice  that  (6.22)  is  an  expression  of  the  following. 

<|>(x)  =  lim  #(x,e)  =  lim  #[x,y-(y-e)] 

£->0  £->0 

,n  >n 


=  lim  ^  Ld2^1L   S   ^(X,y)  . 


This  proves    (6.2.$)   because   if  y  >  a 


,i        i        1    f     (y-a)    f(t)    dt 
t-x  2+y-ar 


—  C3 


The  formula  (6.23)  can  be  expressed  in  a  better  form.  Me   have 


CO 


(6.24)  *(.)  -  i*  A  im±i^ll-n  J/  ^,§1^ 


CO 


f*  o^   1/  ^    n    It   ^  [-i(y-£)1    O"   f     f(t)  dt 

(6.25)  Mx)  =  limQ  -  Im  ^ -i £— -L-  —     t-x-i(y-a)  ' 

Expressed  symbolically,  this  is 

(6.26)  «,)  -  U.  i  JM  .-^^  /  T^SfeS, 

—00 

where  [|  si-  ;  and  y  >  a  .   If  we  take  y  =  a+1  then  (6.26)  becomes 

00 

,r   o^^  U    >    i4    1T   -i(a+l-e)«£  ff(t)  dt 

(6.27)  4>(x)  =  lim  —  Im  e  /  --f  v  r-  • 

£->0  r  1 


■  •'-■•  • 

.  ■ 


■ 


■ 


- 


.' 


. 
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This  means  that  after 


P(x)  .  i  r  4^4fc 


ir  J        t-x-i 

—00 

has  been  found  then  <|>(x)  is  given  by 

<j)(x)  =  lim  Im  P(x-ia-i+€i)  . 

e->0 

As  far  as  the  author  knows,  formula  (6.27)  is  new.   It  appears 
to  be  the  best  of  all  of  the  formulas  we  have  presented  for  the 
solution  of 


f  (X) .  i  r  i±it)dt  ; 

rJ  (t-x)N-a2 


There  is  little  if  any  room  for  an  essential  improvement  of 
(6.27)  as  a  formula  to  be  used  when  f(z)  is  prescribed  analytically 
in  -a  <  Im  2  <  a  ;  and  in  addition  to  this,  the  formula  seems  to 
offer  good  possibilities  for  numerical  analysis  because  it  does  not 
explicitly  contain  a  term  of  the  type  f(z+ia),  as  the  previous 
formulas,  which  involve  one  integration,  do.  Using  observed  values 
of  f(x)  we  can  apply  a  process  of  numerical  integration  to 


CO 


1  /   f(t)  dt 

IT  j  t-X-i 


-co 


and  thus  obtain  P  (x)  as  an  analytic  approximation  to  P(x)  .   The 

a 

values  of  <j)(x)  can  then  be  probed  at  various  depths  a  below  the 
datum  line  by  using 


<|>(x)  =  lim  F  (x-ia-i+is)  . 
e->0  a 

It  should  be  noted  that  [6.23)    can  be  expanded  into 


■ 

'    ■    . 

; 
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2n   >2n 


*(x)  '  ^o  IS  ~wr-  ^  *  J  (x.t)^+(y_a)^ 

+  lim  s=-  [-(y-£).]2?+1  o2n+1,  i  f(y-a)  f(ty  at 

or,  what  is  the  same  thing, 

t6.88)  <,<*) .  urn  f=(-1)?&n)an^BJ  f <y-a!/'t'  a£ 

s->0  feo  l2n)-  Sx2n  ^   (x-t)2+(y-a)d 

+  nm  >sr-  (-i)"!y-^.!2p+1  ag"+1.  .  i  r"(x-t)_  f(t)  at 

In  operational  form  (6.20)  is 

(6.29)   4>(x>  -  lim  cos(y-e)o£>  •  I  /"  (y"a)/(t)  d?j 

a       t        \  £)   ,  x  T  (t-x)  f(t)  dt 
-  lim  sin  (y-e)^  '  =  /  - 'P n    . 

e->0  ^  (t-x)2+(y-a)2 

A  recognition  of  these  expansions  is  helpful  vjhen  we  seek  the 
three-dimensional  analog  of  (6.27). 


■ 


.. 


■ 


. 


- 


• 


■ 


h6 


7.   Polsson's  Integral  Equation  in  Two  Dimensions. 


Let  us  turn  now  to  a  three-dimensional  potential  problem  of 

the  type  discussed  in  Section  2.   If  i!/(x,y,z)    is  a  regular 

potential  function  in  the  half-space  z  >  0;  and  if  ^{x,y ,0+)  =  <}>(x,y) 

where  <{>(x,y)  is  continuous  everywhere,  then  we  know  that  ip(x,y,z) 

is  uniquely  given  by 

00 

(7.1)       V/(x,y,z)  =  i  ff  Li (^ )  d|  d,  ^   ; 

_iJ  t(e,-)2+(Ti-y)2+z2]V2 

If  u.(x,y,a)  =  f(x,y)  is  a  prescribed  analytic  function  then  $>(£,, t)) 
must  satisfy 

.«  [(4-x)  +(n-y)  +a  K'      a  „  0 

which  is  the  counterpart  of  Polsson's  one-dimensional  equation 

(1.1).   As  was  pointed  out  in  the  introduction,  equation  (7«2)  is 

of  considerable  interest  since  it  arises  in  connection  with  certain 

kinds  of  geological  prospecting.  We  proceed  to  give  a  brief 

analysis  of  (7.2)  in  order  to  substantiate  the  fact  that  the 

previous  sections  contain  many  results  which  continue  to  hold  for 

the  two  and  higher  dimensional  generalizations  of  the  equations 

discussed  in  those  sections. 

The  equation  (7.2)  can  be  solved  in  several  different  ways. 

The  method  which  follows  is  based  on  the  semi-group  property  of 

the  operator  0      •   We  define  @ J>   by 

a  a 

03 


/O      M  ■     x      Ff  a  Mf;,n)    dg    dri 

^a  *  a  "Si    J     ri,      T2      '  „XJii/g 


so  that 


[  (e-x)   +(n-y)    +a   ] 


! 


r 

.  •  ) 


M 


.0  ,q  jl  _    i    rr    b rr*  He^)  <%  *\  aex  d^ 

^h^a^  "  1 — 5   /  /  T^ 72 5 — ft  ■*/£•  •'     ?5 3 — 3  T75 

^  _„     [Ux-x)    +(-.:1-y)2+bt]-/2ioo  [(|-C1)2+(-n-rll)24a2]-/2 

=  cn(x,y,b) 

where 

l     rr                   z                         i    rra  ^/^dn  <*i  dTix 
u>(x,y,z)   =  -p-    // 75 ?5 — 5— tt/5 * 5=   //  5 p- — 3  -2/?' 

.»  [Ui-xr-K^-y)2-^2]-/2  2v«>  [(e-y^Ti^f+a2?/2 

The  function  co(x,y,z)  is  a  regular  potential  function  for  z  >  0  and 

.„  [(e-x)s+(n-y)  +a  ]'' 

but  we  also  have 

CO 

,  ftJ,         1     ,,         IP  (z+a)    <J>(e,pJ    d*    dn 

o>(x,y,0+)   =  -p-  lim        /  p  ~g    — ^    ft  '-g/g   i 

27r  z->0-       [(e-x)2+(r)-y)2+(z+a)-]V2 


and   since 


1     /T      (s+a)   4>(€*n)   <%  drt 
27r-i.[(e-x)2+(n-y)2+(ZH-a)2]^2 


is  a  potential  function  in  z  >  0  with  the  same  boundary  value  as 
u>(x,y,z)  we  conclude  from  the  uniqueness  theorem  that 

CO 

a)(x,y,z)   =   p-    //  S — P      '^'a' — 5 — ^   */2    • 

J_co   [(^-x)2+(^-y)2+(z+a)2]-/2 

Therefore 

CO 

,0    /O     1          1          ui        1      rr           (b+a)    <MS»")   d|   dn 
^h^a   "P  =  «"(x,y,b)  =  -p-   .//  \g  ,  T  .6  ! —  »a"*/g 


or 


(7-3)  &,&♦-&*♦■ 


- 


• 


' 


I 


• 


d 


■ 


■ 

- 


( 


• 


' 


■ 
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The  property  (7.3)  can  be  used  to  solve  (7.2).   .'ith  0  <  z  <  a, 

n 

the  application  of  (y     to  each  side  of  (7.2)  produces 

00 

00  CO 


=  •&  ff  fv  e~(a+z)r  j    (rp)   dr  4>(C,r:)    d^   dn 


2ir  . 

-00  0 


2 g 

where  p  =  y(£-x)  +(rj-y)   .   The  last  equations  is  the  same  as 

03  CO 

-co"      0 


-  —    /T  (a-z)    i(^n)    d^    dn 

27f  J_i   [U-x^^-y^a-z)2]^ 


If  we  let   z   ->  a-0  we   find 

CO  CO 


<j>(x,y)   =   lim       -^JJ  j    r   e"ar.    2  X"  Tfe~Tr     J0(rp )dr   M^MI*! 

—00     0 

-0.  * 

Since 

7%<rp)  .   (^  +a,VT0(rp)  -  [-*  +i^]J0(rp) 
ox         oy  op 


a 


,2,     ,„     N    _     ,d2  S2     v,     ,        ,  r     d2 


=   -  r2  JQ(rp) 
it   follows   that 


V2n  j  (ro)   =    (-Dn  r2n  JQ(rp) 


and   hence 

CO  CO 


4>(x,y)   =   lim       2^I^4nTT     ^  h  ff  [  re'^¥r^    dr  K.rifc  *i 
z—>,a-"        n— 0  /-\ 


.  : 


•■ 


: 


. 


. 


. 


•  t 


- 


)  ■ 


■: 


«>      f    ,,n   2n_2n      ^ 
4>(x,y)  .  iia       2JZ   ('11/ti,        dX   k  -   tZ 


z->a-        n=0 
which  yields 
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(7.M  <|>(x,y)   -  lim       2  ^Z       |Lf«     V        f(x'y) 

z->  a-        n=0        * 


00 


1     pf        a   ff^ri)    d£    dvi 

"  s-  //  [(e.x)^(,:y)^]'/-- 

for  the  solution  of  (7-2). 

We  can  now  ask  for  verions  of  (7.4)  which  can  be  used  to 
compute  cj)(x,y)  when  f(x,y)  is  a  function  derived  from  experimental 
observations.   For  numerical  purposes  there  are  many  approximations 
that  can  be  used  in  place  of  (7.4)  and  these  approximations  can  be 
deduced  by  using  the  methods  of  Section  6.   In  the  main  part  of 
what  follows  we  will  confine  ourselves  to  the  derivation  of  a 
formula  which  can  be  based  on  an  iteration  procedure  analogous 
to  that  presented  in  the  last  section. 

The  solution  of  (7.2)  can  be  expressed  in  the  symbolic  form 


*=   1   f  = 


'a     l-(l-<9a) 

and  from  this  we  find  the  formal  symbolic  solution 


If  we  suppose  that 


tee    a 


n  K 


(7.5)  *n  =  g  a-^r  f 

is  known  then  <j>  ,  is  given  by  the  recurrence  formula 

i    rr      a  bnli,r\)  d£  at) 

<L.-,(x)  =  f(x)  +  <b  (x)  -  4-  //  i 2     2  V2  * 

n+1  n       2   Jj     [(e-x)2+(n-yr+aZ]'/2 


■ 


■ 
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This  also  means  that  <f>  must  satisfy 


n 


(7.6)  4=  //  § 2  g  ^./g  =  f  -  (4>n+1-4>n) 

-«  [(e-x)2+(-,;-y)2+aV/2         n+1  n 

=  f  -  (!-<?  )n+1  f . 

cL 

As  in  Section  6  it  can  be  shown  that  if,  for  example,  <{)  =  '-/5<f>-i  with 

<}>,  absolutely  integrable,  then  i  converges  uniformly  to  <\>   as  n->«  , 

Use  of  the  semi-group  property  of  (/     allows  us  to  write  (7.6) 

a 

as 

C9a+n=((n+l)tfa-i£i|^^a...  <-Dn<?(n+1)a}f 

from  which 

(7.7)  in(x,y)  =  (n+l)f(x,y)  -  i^i2  ^&  f  ...  (-l)n^a  f  . 

If  we  use  polar  coordinates 

£    -   x  =  p  cos  9  r\   -  y  =  p  sin  0 

f(^*n)  =  f(x+o  cos  9,  y+p  sin  9)  =  F(p,©;x,y) 

then  (7.7)  becomes 


(7.8)   4>Jx,y)  =  (n+1)  f(x,y) 

(n+l)n     a     .  (n+l)n(r>-D.    2a 


n 

002.  \    --;-^r-  (py,*/^        3J      (^7^/ 


0  o  (_1)  n£l 


(.       •  •  •  (77nW72  y' 


This  can  be  expressed  in  essentially  the  same  form  as  its 
one-dimensional  counterpart  if  we  notice  that 


•  - 


.  ) 


: 


- 


. 


■ 


. 


■ 


■• 


• 


• 


■ 


• 


• 


■ 


. 


• 
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GO 


— * — 0  *   -,,*  =    /    r   e  j    ( rp )    dr 

[p2+K2a2]5/2      ^ 


00  00 

t   #™t        2    r  s^  rx  dr       2  Tm    f      eirx  dt 


0  ^  j         ^_ 7T  J       ^^ 2 

pTT-p  P7/T-P 


00 


*-* ^Re    f  — i~  4    /"  e-'-^e1"  dr  dT 


J  p  n  -p  o 


Hence 


<!>n(x,y)  -    (n+1)    f(x,y) 


oo    2tT  oo  oooo 

1     f  r    _  ._    .         _      f       1  c)      P    irTf-(n+l)n     -ar      *,».»*» 

--5  jjp  F  dO  dp    •  Rej——^!   e       {  — gj—  e      ...}drdT   . 

T  o  o  p  n/x  -pd       o 


The   integral 


(n+l)n     -ar        (n+l)n(n-l)      -2ar 
— oT —   e  H  =T! ' 


f  eirx{  ...  )    dr 


_.  ,   .,n     -nar 

0  .  . .    (-1)      e 

was  evaluated  in  Section  6  and  it  is  equal  to 

Tm   no  f     1   .  (i)n  (n+1)  I    i 

Im  =  Re  [  -  ^^   +  a(p.i)p(p+i)...(p+ni)  3 

■ 

where  £  =  x/a.  We  therefore  have 

oo    2tt  oo 

(7.9)     <L(x,y)   =    (n+1)    f(x,y)   +i   f(PFd0do    [      1        £.(-*£& 

*     0  J0  p  7^         a  ** 

oo    2tt  °° 

(n+D  I  r  r~  „  ~  „„  ..  r i a 


2 


p   F  d0  dp  Re       = £=0^(P)    dr 


a7r"      (T  cT  p'  Yt"2^ 


where 


.    " 


. 


;      •       • 


. 


• 


■ 
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4i«\       (i)"  . 

°  {^'    '    (p-i){3(  ,4-i)T3+2i)  .  .  .  O+ni) 


The  form  of  equation  (7.9)  can  be  improved  by  noting  that 

CI 

O 

d     i      a 


1  o     /      a *     ,  7T  a 

y^2  6'      a2+r2  2   [p^a2]^2 


and  that 


oo    27T  oo 

a7r        0  0  p  yr-Pd 

=      -    (n+1)    +  2    . 


The  use  of  these  results  allows  us  to  write 

oo  2ir 
(7.10)     *n(x,y)  =  2  f(x,y)   -    *.  f  f      »  gp  ffif 

0  0       Lp  +a  J 


co    2tt 


_ln+l)!     r  [  [P.f(x,y)]0    d9   dP   Re    f        *  4   ^(j3)    dx 

<-  /  s  o  'J      ^  *■ 

a?r       ~  ^  ^/-^   _c 


0   0 


^— ^  on 
p    1/t  -o 


or 


(7.11)      4> 


in  e)p   dp    d9 


2tt   oo 

/        \        o^/        <>         1     /~   -r  a   f  (x+p   cos  ©,y+o      s 
n(x,y)   =   2   f(x,y)    -  ^  J       L_C >J 

0  0  [p     d   J 

2ir   oo 
^ n"!o ^  l    //    [f(x+ao   cos  9,y+ap   sin  Q)-f(x,y)]o   dp   d©' 

0      0  03 


•  ^~~ 2   °£ 

?    70  -P 


which  may  be  useful  for  numerical  analysis  if 


p  Vp2-p2 


is   tabulated 


) 


■ 


.      ■  ■      1 


I 


' 


: 


5? 

Let  us  conclude  with  a  derivation  of  an  analog  of  (6.27). 
The  unknown  function  <f>(x,y)  can  be  expressed  as 

<|>(x,y)  =  lirn  ^(x,y,e)  =  lim  ^[x,y,z-(z-e)] 
e->0  e->0 


00 


I-(z-e)]n   an 

e->0  n=0 

which  implies 


-   llm^  ^__  sn — =     #(x,y,z) 

e->0  n=0  n*  dz 


(|>(x,y)   =   lim    ^  H2-e 

e->0  n=0  n 


1^   j£.  1    ;T       (z-a)  f^^)  gg  dn 

CO 

=   lim     e"(2"e>^.  J-    /T         (z-a)    f(g.tin)    dg    dy 
e ->0  27f  ^  [(S-x)*+(T)-y)2+(z-a)2]'/2 

where  ^  =  -^_  and   z  >   a    .      This   can  be  expanded   into 

e->0  fen        <2n>!   S^       »  W    [(e-xl^fn-yl^fz-a)^'^ 

.  lim  2Z  ("1,)^(!Te)2n  (-4  +  -&ni  r"      (z-a)  fti,aj  dg  an 
£_>o  fen  — rssn       (^    sp>'   sf  j_;  [u.x)^+(n.J^;(2.a).).v^ 

+  lim  f~   (-D^tz-e)5"1-1  ,_&!  +  _df.n+l  1  ff      fje.i)  g  d. 


,n,  ,,   .2n   ^2    ^2 


In  particular  if  we  take  z  =  a+1 

♦(x.y)  .  iim  ■£-  (-1)n[^T£)gn  (^  +  4>n^  ff   ^  f  %g 

e->On=C     l2n)-     3X2   Sy2   ^'M-xWn-yl^l]'72 
til.  f  (-l)n+1(a+l-£)an+1  ,  3s    Ss.n+1  1  /T  fU.n)  «  dn 


: 


■ 


■ 


s 


-   ' 


• 


- 
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- 
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■ 

r 
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■ 


■ 
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